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The Input-output characteristics of an ultrasonic testing system used for 
"stress wave factor" measurements were studied by coupling the transmitting 
and receiving transducers face to face, without a specimen In between. Some 
of the fundamentals of digital signal processing used are summarized. 

The Inputs and outputs were digitized and processed In a microcomputer by 
using digital signal-processing techniques. The entire ultrasonic test system, 
Including transducers and all electronic components, was modeled as a discrete- 
time, linear, shift-invariant system. A digital bandpass filter was Introduced 
to reduce noise effects on the output signal. The output due to a broadband 
Input was deconvolved with the Input to obtain the unit-sample response and the 
frequency response of the discrete-time system. Then the Impulse response and 
frequency response of the continuous-time ultrasonic test system were estimated 
by Interpolating the defining points In the unit-sample response and frequency 
response of the discrete-time system. The ultrasonic test system was found to 
behave as a linear-phase bandpass filter. 

The unit-sample response and frequency response of the discrete-time model 
of the test system were used to compute the output of the test system for a 
variety of Inputs. The agreement between predicted and measured outputs was 
excellent for rectangular pulse Inputs of various amplitudes and durations and 
for tone burst and single-cycle Inputs with center frequencies within the pass- 
band of the test system. The Input-output limits on the linearity of the 
system were determined. 

INTRODUCTION 


Conventional ultrasonic testing (UT) Is conducted either In the through- 
thickness transmission or pulse-echo modes (ref. 1). In the through-thickness 
transmission mode, the transmitting and receiving transducers are coupled to 
opposite faces of the structure under Inspection, and the transmitted wave 
field Is analyzed. In the pulse-echo mode, the transmitting and receiving 
transducers, which may be combined Into a single transducer, are coupled to the 
same face of the structure under Inspection, and the reflected wave field Is 
analyzed . 

Recently Vary et al. (refs. 2 and 3) Introduced an ultrasonic nondestruc- 
tive evaluation (NDE) parameter called the stress wave factor (SWF). It Is 
similar to the pulse-echo test mode In that separate transmitting and receiving 
transducers are coupled to the same face of the structure. However, unlike 
conventional pulse-echo testing, which Is generally limited to the analysis of 
nonoverlapping reflected wave echoes, the SWF Is also valid for the analysis 
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of overlapping echoes. Specifically, an Input pulse having a broadband fre- 
quency spectrum Is applied to the transmitting transducer and the number of 
oscillations of the output signal at the receiving transducer exceeding a pre- 
selected voltage threshold Is defined as the SWF. The SWF has been correlated 
with mechanical properties of carbon-fiber-reinforced composites (refs. 2 
to 4). Williams et al . (ref. 5) theoretically and experimentally studied the 
ultrasonic Input-output characteristics of the SWF test configuration In a 
thick, Isotropic, elastic plate. The extension of that work to thin plates Is 
currently In progress (ref. 6). 

An Important step toward the quantitative analysis of Input-output rela- 
tions In any ultrasonic NDE procedure Is the quantitative characterization of 
the experimental test system without a test specimen. The effects on the out- 
put signal of the ultrasonic transducers, the coupling of the transducers to 
the test specimen, and electronic components such as filters, amplifiers, 
attenuators, and cables In the experimental system must be characterized before 
the test specimen's effects on the output signal can be Isolated. 

This study Is part of an overall effort to develop quantitative analyses 
of the SWF and computer-aided nondestructive evaluation (CANDE, pronounced 
"candy") capabilities. The experimental UT system Is characterized by directly 
coupling the transmitting and receiving transducers face to face without a test 
specimen. Input and output signals are digitized with a digital oscilloscope 
and are processed with a microcomputer by using digital signal-processing tech- 
niques. The transfer function of the experimental UT system without any test 
specimen was obtained. 

The results of this study should provide a useful example In the charac- 
terization of any UT system that has separate transmitting and receiving trans- 
ducers. Furthermore, developments In CANDE should be facilitated by the 
digital signal-processing procedures summarized herein. 


FUNDAMENTALS OF DIGITAL SIGNAL PROCESSING 

A few results In digital signal processing are summarized In this section. 
Digital signal-processing techniques are discussed extensively In the litera- 
ture (refs. 7 to 10). Primarily, the notations In reference 7 are followed In 
this outline. 

Digital signal processing Is concerned with the representation of signals 
by sequences of numbers and the processing of those sequences. Sequences 
correspond to discrete-time signals derived from sampling continuous-time 
signals. The notation x(n) denotes a sequence of numbers whose entries 
depend on the Independent parameter n. The sequence x(n) Is defined only 
for Integer values of n and represents successive samples of a continuous- 
time signal. For example, the sequence x(n) Is derived from periodic 
sampling of the continuous-time analog signal x a (t), where t Is time, 
according to 


x(n) = x a (nT) (1) 

where T Is the sampling period. The reciprocal of T Is the sampling rate 
or sampling frequency. The availability of high-speed digital computers and 
efficient signal-processing algorithms has accelerated the Implementation of 
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digital signal-processing techniques over that of continuous-time, analog 
signal-processing techniques. With the proper analysis, results from the 
digital signal processing of a sequence derived from sampling a continuous-time 
signal can accurately approximate results from continuous-time signal analyses. 

A discrete-time system Is defined mathematically as a unique transforma- 
tion that maps an Input sequence x(n) Into an output sequence y(n). 
Discrete-time, linear, shift-invariant systems are discrete-time analogs of 
continuous-time, linear, time-invariant systems. Table 1 summarizes the Input- 
output properties of discrete-time, linear and linear shift-invariant systems 
and their corresponding analogs In continuous-time systems. Discrete-time 
systems can be further characterized as stable and causal. A stable system Is 
one for which every bounded Input x(n) produces a bounded output y(n). A 
causal system Is one for which the output y(n) for n equal to ng depends 
only on the Input x(n) for n less than or equal to ng. All subsequent 
discussions here are limited to discrete-time, linear, shift-invariant systems 
that are stable and causal. 

The response of a linear, shift-invariant system can be characterized by 
the unit-sample response h(n) via the convolution sum given In equation (1-1) 
In table 1. As a consequence, It can also be shown that the steady-state 
response of a linear, shift-invariant system to a sinusoidal Input Is sinusoi- 
dal of the same frequency as the Input but with a magnitude and phase deter- 
mined by the system. It Is this property of linear, shift-invariant systems 
that makes representations of signals In terms of sinusoids or complex exponen- 
tials (i.e., Fourier representations) so useful In linear system theory. Also, 
It should be mentioned that the process of obtaining y(n) from known x(n) 
and h(n) is called convolution, and the process of obtaining h(n) from 
known x(n) and y(n) Is called deconvolution. 

Table 2 summarizes Fourier representations of continuous-time signals and 
of sequences. Fourier representations appear as equation pairs whose constitu- 
ents are often referred to as the synthesis and analysis equations, as Indi- 
cated In table 2. In equation (2-lb), X a (j«) is a continuous function In 
radian frequency Q and Is the Fourier transform of the continuous-time signal 
x a (t), where the symbol j Is defined as the square root of -1. If x a (t) 
has units of voltage, X a (jG) has units of voltage per radian frequency (or 
simply volt-second); X a (j«) Is In general complex and can be specified by Its 
magnitude and phase. 

If x a (t) is a real-valued function of time, the magnitude and phase of 
X a (jfi) will be an even and odd function of frequency, respectively. Thus, it 
Is only necessary to present the magnitude and phase of X a (jQ) for positive 
frequencies. If x a (t) Is real and an even function of time, the phase of 
X a ( jn) Is zero. If x a (t) Is real and an odd function of time, the phase 
of X a (jn) Is ±ir/2. If x a (t) is shifted (I.e., advanced or delayed) In 
time, the magnitude of X a (jft) Is not affected, but a phase linearly propor- 
tional to frequency (i.e., linear phase) Is Introduced In the phase of X a (ja). 
More specifically, the phase (In radians) Is proportional to the radian fre- 
quency with a proportionality constant having units of time. The negative of 
the value of the proportionality constant Is the time delay (ref. 9). 
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The remainder of table 2 deals with the Fourier representations of 
sequences. If properly computed, the Fourier transform X a (jQ) can be 
approximated by the Fourier representation of a sequence derived from sampling 
x a (t). Examples of the Fourier representations of continuous-time signals 
and sequences are Illustrated In figure 1 and will be discussed shortly. 

In equation (2-2b) In table 2 X(e3 w ) Is a continuous function of radian 
frequency u and Is the (discrete time) Fourier transform of a sequence x(n). 
The frequency o> Is In units of radians per Increment of n, which Is simply 
radians. If x(n) has units of voltage, X(e3 w ) has units of voltage per 
radian. The Fourier transform of a sequence Is useful for analyzing general 
sequences. The discrete Fourier series representations of a perlodlc^sequence 
x(n) of period N are also shown In table 2. A periodic sequence x(n) of 
period N Is one such that x^n) equals ?(n + N). In equation (2-3b) X(k) 
Is a periodic sequence of period N, and the numbers In this sequence are the 
discrete Fourier series coefficients of the periodic sequence x'(n) having 
period N. If x(n) has units of voltage, ^(k) has units of voltage also. 

The discrete Fourier series of a periodic sequence serves as a prelude to the 
discrete Fourier transform (OFT). In equation (2-4b) X(k) Is a sequence of 
length N and Is called the DFT of a finite-length sequence x(n) of length 
N (ref. 7). If x(n) has units of voltage, X(k) has units of voltage also. 
The OFT Is useful for digital signal processing because sequences processed by 
a computer are of finite length and also because efficient OFT computational 
algorithms are available. 

By comparing the entries In table 2 for the discrete Fourier series of a 
periodic sequence of period N and the OFT of a finite-length sequence of 
length N, It Is observed that the DFT representation Is obtained from the 
discrete Fourier series representation by Interpreting the finite-length 
sequence as one period of a periodic sequence. The properties of the OFT are 
similar to those of the Fourier transform, except that because of the Implied 
periodicity, shifts of x(n) In n by one period N and an Integer multiple 
of one period are Indistinguishable and a shift In n of larger than N Is 
the same as a shorter shift. 

The DFT was computed via a fast Fourier transform (FFT) algorithm. The 
Inverse DFT computed via an Inverse FFT (IFFT) algorithm was used to evaluate 
x(n) from X(k). The FFT computation Is particularly efficient when the 
length of the sequence Is an Integer power of 2. A sequence of length N Is 
called an N-poInt sequence. The DFT of an N-poInt sequence Is called an 
N-poInt DFT and Is computed via an N-poInt FFT algorithm. Figure 1 Illustrates 
that the Fourier transform of a continuous-time signal can be reconstructed 
approximately from the DFT of a sequence derived from sampling the continuous- 
time signal. Again, note that where DFT relations are concerned, a finite- 
length sequence Is represented as one period of a periodic sequence. 

Figure 1 shows Fourier representations of continuous-time signals and of 
sequences derived from sampling the continuous-time signals. In part (1) of 
figure 1(a), the continuous-time signal x a (t) Is assumed to be bandllmlted 
In frequency with Fourier representation X a (jft) as shown. The highest fre- 
quency of x a (t) Is assumed to be Si 0 /2. In part (11) of figure 1(a), a 
sequence x(n) results from the periodic sampling of x a (t) with a sampling 
Ing period of T. The Fourier transform of the sequence x(n) Into a contin- 
uous function X(e3«) Is also shown. It Is observed that X(eJ“) can be 
obtained from X a (;)G) by the superposition of an Infinite number of X a (jfl) 
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shifted In frequency u by Integer multiples of 2ir. The magnitude and fre- 
quency of X(e^ w ) are scaled from those of X a (j«) by 1/T and T, respec- 
tively. The Fourier transform X a (jQ) can be recovered exactly, except for 
the scaling factors, from the Fourier transform X(eJ“) of the sequence x(n) 
In the Interval of -ir < w < ir by low-pass filtering X(eJ u ) with a cutoff 
frequency of ir, on one condition. The condition Is that Q 0 T/2 < ir so that 
there will be no overlapping of superposed X a (jQ) In forming X(eJ w ). This 
condition Is known as the "sampling theorem" and assures that If a continuous- 
time signal x a (t) Is sampled at a frequency greater than twice the highest 
frequency of x a (jQ), then X(eJ“) Is Identical to X a (jn) In the Interval 
-ir < w < ir, except for scaling factors. Because X a (jQ) Is recovered, the 
continuous-time signal x a (t) can be recovered from the sequence x(n). This 
minimum required sampling frequency Is called the Nyqulst rate or the Nyqulst 
frequency (refs. 7, 9, and 10). The distortion due to overlapping of super- 
posed X a (jfi) when the sampling theorem Is violated Is called aliasing. In 
part (11) In figure 1(a), the case where there Is no overlapping of superposed 
X a ( jn) Is Illustrated; thus the signal Is not aliased. 

Whereas figure 1(a) deals with general continuous-time signals, 
figure 1(b) deals with finite-duration continuous-time signals. In part (1) 

In figure 1(b), a continuous-time signal x a (t) with nonzero values over a 
finite time duration and Its Fourier representation X a (jQ) are shown. The 
duration where the signal Is nonzero Is assumed to be t 0 , and sampling the 
continuous-time signal with a sampling period T Is assumed to result In a 
sequence of length N. It Is observed that X a (jQ) for a finite-duration 
signal Is not bandllmlted In frequency. In part (11) In figure 1(b), a 
sequence x(n) results from the periodic sampling of x a (t) with a sampling 
period of T. All numbers In the sequence x(n) are zero for n < 0 and 
for n > N - 1. Thus x(n) Is a finite-length sequence of length N. The 
Fourier transform of the sequence x(n) Into the continuous function X(eJ w ) 

Is also shown. Similar to the results In figure 1(a), X(ei“) can be obtained 
from X a (jQ) by the superposition of an Infinite number of X a (jn) shifted In 
frequency u by Integer multiples of 2ir. However, for this case, because 
X a (jfi) Is not bandllmlted In frequency, there Is overlapping of superposed 
X a (jfl) In forming X(e^ w ). Thus there Is aliasing and X a (jft) cannot be 
recovered exactly from X(eJ w ) In the Interval -*<«<*. However, If the 
sampling frequency Is greater than twice most of the significant frequency 
components In X a (ja), then X a (jn) can still be approximated adequately by 
X(e3 w ). The DFT X(k) of the finite-length sequence x(n) Is often computed 
Instead of X(eJ w ). However, In DFT evaluations, a finite-length sequence 
Is represented as one period of a periodic sequence. 

In part (111) of figure 1(b), a periodic sequence 'x(n) of period N 
constructed by using the N-poInt sequence^ 'x(n) as a period Is shown. The 
discrete Fourier, series representation X(k) of the periodic sequence x'(n) 

Is also shown; )f(k) Is a periodic sequence of period N. The DFT X£k) of 
the finite-length sequence x(n) can be Interpreted as a period of )T(k) for 
jO<k<N-l as shown. It Is observed that X^ 01 ) can be obtained from 
X(k) by Interpolating points In ”)f(k). Because X a (jQ) can be approximated 
adequately by X(e^ w ) If the sampling frequency Is greater than twice most of 
the significant frequency components In X a (jn), Interpolations of X(k) or 
X(k) can be used to estimate X a (jQ). The numbers In X(k) are scaled 
approximately from the values of X a (jtt) by a factor of 1/T, and the spacing 
between successive numbers In the sequence X(k) can be Interpreted as 2-ir/( NT) 
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In radian frequency, which Is the frequency domain resolution. Again, the 
reconstruction of the Fourier transform of a finite-duration, continuous-time 
signal from the DFT of a finite-length sequence by Interpolation Is generally 
an approximation due to aliasing. 

Table 3 lists the conditions for the special case where there Is exact 
correspondence (to within a scaling factor) of the Fourier transform of a 
continuous-time signal and the DFT of a finite-length sequence derived from 
sampling (ref. 10). Also, the N-poInt sequence x(n) can be considered as an 
M-poInt sequence If M Is greater than N with the last M-N numbers In the 
M-poInt sequence being zero. As discussed here next regarding convolution 
procedures, sometimes It Is convenient to extend the length of a sequence by 
appending zeroes. 

Because the response of a linear, shift-invariant system can be evaluated 
via the convolution sum In equation (1-1) In table 1, convolution representa- 
tions are discussed here. Table 4 summarizes convolution representations of 
continuous-time signals and of sequences. The Fourier representations of the 
convolutions are also given In table 4. Equation (4~la) gives the linear con- 
volution of continuous-time signals x a -|(t) and x a 2 (t), resulting In x a 3 (t). 
The Fourier representation of the convolution Is given In equation (4-lb). 
Convolution In the time domain results In multiplication In the frequency 
domain. Conversely, because of the duality between time and frequency domains 
In Fourier representations, convolution In the frequency domain results In 
multiplication In the time domain, such as time windowing (ref. 9). The con- 
cepts of periodic and circular convolutions are Introduced as preludes to the 
efficient evaluation of linear convolution of sequences using DFT procedures. 

Equation (4-2a) In table 4 gives the periodic convolution of two periodic 
sequences ^(n) and x^n), each of period N, resulting In a periodic 

sequence X 3 (n) of period N. By restricting attention to one period of the 
periodic convolution, circular convolution of two N-poInt sequences x-j(n) 
and X 2 (n) to obtain an N-poInt sequence X 3 (n) can be defined from the 
periodic convolution and Is given In equation (4-3a). The two periodic 
sequences x^n) and Y 2 ( n )’ each of P er1od N > are constructed by using the 

N-poInt sequences x-|(n) and X 2 (n) as periods, respectively. Equations (2-3b) 
and (2-4b) In table 2 have been used to obtain the simplified equation (4-3b) 
because ^(n) and x^) are equal to x-|(n) and X 2 (n) for n ranging 
from zero to N - 1. Thus X^k) and X 2 (k) for 0 < k < N - 1 become X-|(k) 
and X 2 (k), respectively. The circular convolution of N-poInt sequences Is 
called an N-poInt circular convolution. The designation "circular" Is derived 
from a graphical representation of how an N-poInt sequence can be used to con- 
struct a periodic sequence of period N (ref. 7). It can be Imagined that the 
N points from the N-poInt sequence are equally spaced In angle around a circle 
with a circumference of exactly N points. The periodic sequence of period 
N Is obtained by traveling around the circumference of the circle a number at 
a time recording the N-poInt sequence repeatedly. Also, a rotation of the 
circle corresponds to circular shifting of the N-poInt sequence (ref. 7). 

The linear convolution of an N-|-po1nt sequence x-|(n) and an N2-po1nt 
sequence X 2 (n) Is defined as 


x 3 (n) 


X) x,(m) x 0 (n - m) 
m=0 1 1 


( 2 ) 
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Because the resulting sequence X 3 (n) Is of length N] + N 2 - 1 (ref. 7), 
the evaluation of linear convolution via circular convolution requires all 
sequences to be of the same length, which has to be N] + N 2 - 1 . 

Equation ( 4-4a ) In table 4 gives the linear convolution of finite-length 
sequences x-|(n) and X 2 (n) to obtain a finite-length sequence X 3 (n) 
evaluated via circular convolution. Computation of circular convolution using 
DFT Is preferred over the direct evaluation of linear convolution because of 
computational efficiencies. The (N-| + N 2 - l)-po1nt DFT's of the two convolv- 
ing sequences are computed via an FFT algorithm and then multiplied according 
to equation ( 4- 4b ) . The (N-) + N 2 - l)-po1nt Inverse DFT of the product Is 
computed via an IFFT algorithm, and the result Is the desired linear convolu- 
tion. Because the original sequences x-|(n) and X 2 (n) are of lengths N-| 
and N 2 , respectively, N 2 - 1 and N-| - 1 zeroes are appended to x-](n) 
and X 2 (n), respectively, to yield two (N-| + N 2 - l)-po1nt sequences for the 
(N-| t N 2 - 1 )-po1nt circular convolution. Alternatively, If linear convolu- 
tion of two ( N] + N 2 - l)-po1nt sequences Is conducted via an (N] 4- N 2 - 1 ) - 
point circular convolution, the last N 2 - 1 numbers In one sequence to be 
convolved must be zeroes and the last N-j — 1 numbers In the other sequence 
to be convolved must be zeroes. In practice, the sequences are extended to a 
length that Is greater than Ni 4 - N 2 - 1 and Is an Integer power of 2 for 
efficient FFT computation. 

Figure 2 shows convolution representation of continuous-time signals and 
their sequences derived by sampling. In figure 2(a), continuous-time signals 
x a i(t) of duration t-| and x a 2 (t) of duration t 2 are convolved to form 
the continuous-time signal x a 3 (t) of duration t] 4- t 2 - In part (1) of 
figure 2(b), the continuous-time signal x a ](t) Is sampled with a sampling 
period T and Is assumed to produce an N-j-poInt sequence. Then N 2 - 1 
zeroes are appended to form an (N] 4- N 2 - l)-po1nt sequence x-|(n). Next a 
periodic sequence ’x'i(n) of period N-] + N 2 - 1 Is constructed by using 
xi(n) as a period. Similarly, In part (11) In figure 2(b), the periodic 
sequence X 2 (n) of period N] + N 2 - 1 Is obtained from x a 2 (t), which 
Is assumed to produce an N 2 -po 1 nt sequence when sampled with a sampling period 
T. Periodic convolution of x'i(n) and ^(n) results In the periodic 
sequence ^(n). Then x 3 (n) can be Identified from the circular convolution 
portion of x^n) by restricting attention to one period of the periodic 
sequence. It Is observed that x a 3 (t) can be approximated from the Inter- 
polation of the (N-| 4- N 2 - l)-po1nt sequence x 3 (n). The numbers In X 3 (n) 
are scaled approximately from the value of x a 3 (t) by a factor of 1/T and the 
spacing between successive numbers In the sequence X 3 (n) can be Interpreted 
as T. Because the circular convolution Is obtained via a DFT, restrictions on 
the DFT apply to the convolution also. Thus the sampling theorem applies to 
sampling the continuous-time signals If x a 3 (t) Is to be approximated 
adequately by the sequence x 3 (n). 

In practice. If a continuous-time signal Is to be sampled, the sampling 
frequency and the sampling length must be specified. Table 5 summarizes the 
significance of the sampling parameters. Equation (5-1) Is based on the 
earlier discussion In relation to figure 1(a). 


EXPERIMENTAL SYSTEM 

A schematic of the experimental system Including the specimen Is shown In 
figure 3. This Is a typical stress-wave-factor test configuration. The goal 
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of this study was to characterize the experimental test system without the test 
specimen. Thus, all testing In this study was conducted with the transmitting 
and receiving transducers directly coupled face to face without any test 
specimen In between. 

The system consisted of a pul se/f unction generator (Wavetek model 191); a 
5-MHz (corresponding to -3-dB point) low-pass filter (Allen Avionics F2516); 
broadband (0.1 to 3.0 MHz) transmitting and receiving transducers (Parametrlcs 
VI05) having an approximately flat sensitivity of -85 dB relative to 1 V/pbar; 
an ultrasonic Interface couplant (Acoustic Emission Technology SC-6); an ultra- 
sonic preamplifier (custom built by Parametrlcs) having a gain of 40 dB In the 
frequency bandwidth 1 kHz to 10 MHz; a variable-frequency filter (A.P. Circuit 
Corporation AP 220-5), which could be used as a low-pass, high-pass, bandpass, 
or bandstop filter In the frequency range 10 Hz to 2.5 MHz; a digital oscillo- 
scope (Nlcolet model 2090 with plug-in model 204-A), which could sample and 
store analog signals at sampling frequencies from 0.05 Hz to 20 MHz (corre- 
sponding to sampling periods from 20 sec to 50 nsec) up to 4096 points; and 
an IBM personal computer (IBM PC), which was Interfaced with the digital 
oscilloscope. 

The variable-frequency filter was used as a 0.4- to 2.6-MHz (corresponding 
to -3-dB points) bandpass filter, selected on the basis of the receiving trans- 
ducer frequency response. All signals were sampled at a sampling frequency of 
20 MHz (l.e., one sample every 50 nsec) for 256 points. Thus each recorded 
signal was 12.8 ysec In duration. The 5-MHz low-pass filter was selected via 
antialiasing considerations. Table 6 Is a BASIC language program listing for 
the FFT and IFFT computations used In this study. All DFT and Inverse DFT 
evaluations In this study were conducted by using this program. 


DIGITAL CHARACTERIZATION OF EXPERIMENTAL SYSTEM 

The objective of this study was to obtain the system transfer character- 
istics between the Input to the transmitting transducer and the output from 
the receiving transducer after signal conditioning. A schematic of the experi- 
mental system Is shown In figure 3. Once again, note that the experiment was 
conducted with the transmitting and receiving transducers coupled face to face 
without a test specimen. Thus, the "system" under study consisted of the 
transmitting transducer, the coupling between transmitting and receiving trans- 
ducers, the receiving transducer, the ultrasonic preamplifier, the variable- 
frequency filter, and all necessary cables. 

A continuous-time system can be described by Its Impulse response h a (t) 
and Its frequency response H a (jQ). In this study, the h a (t) and H a (ja) 
of the test system were computed from the unit-sample response h(n) and the 
DFT of h(n) denoted H(k), respectively, of a discrete-time model of the 
continuous-time system. 

The experimental system was assumed to be a discrete-time, linear, shift- 
invariant system. As shown In figure 3, channel 1 of the digital oscilloscope 
sampled the continuous-time Input signal to the transmitting transducer. The 
resulting discrete-time signal Is called the Input sequence and Is denoted by 
x(n) . Similarly, channel 2 of the digital oscilloscope sampled the continuous- 
time output signal from the receiving transducer after preampllf Icatlon and 
bandpass filtering. The resulting discrete-time signal Is called the output 
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sequence and Is denoted by y(n). Thus, the objective of this study was to 
obtain the discrete-time system transfer characteristics between x(n) and 
y(n). The discrete-time system besides being linear and shift-invariant was 
observed to be stable and causal. A schematic of the digital signal-processing 
procedures applied to obtain the unit-sample response and frequency response 
of the discrete-time system Is shown In figure 4. 

The pulse/functlon generator produced an approximately rectangular pulse 
20 nsec In duration. The rectangular pulse was passed through the 5-MHz low- 
pass filter before It reached the digital oscilloscope. Figure 5 shows the 
time histories of the Input and output signals by linearly Interpolating the 
defining points In the Input and output sequences x(n) and y(n). The Input 
and output signals were approximately 3 and 4 ysec long, respectively. Thus 
for the sampling period 50 nsec, only approximately the first 60 and 80 numbers 
In the sequences x(n) and y(n), respectively, were nonzero. Actually, the 
sequences x(n) and y(n) were of length 256 after sampling; that Is, the 
digital oscilloscope recorded 256 points each for the Input and output signals. 
Then 768 zeroes were appended to the sequences to produce x(n) and y(n) of 
length 1024. Extended sequences were used to facilitate evaluation of linear 
convolution via the DFT, as discussed later. The DFT's X(k) and Y(k) for 
x(n) and y(n), respectively, were evaluated by using a 1024 -point FFT algo- 
rithm based on 1024-point x(n) and y(n). Because the Input and output 
signals were of finite duration, they were not bandllmlted In frequency. How- 
ever, because the sampling frequency (20 MHz) was greater than most of the 
significant frequency components In the Input signal (significant up to 5 MHz) 
and the output signal (significant up to 2.6 MHz), the sampling theorem assured 
an adequate approximation of the Fourier transform of the Input and output 
signals by the DFT of the Input and output sequences. Approximations of the 
magnitude and phase of the Fourier transform of the Input and output signals 
are shown In figure 5. The magnitude and phase of the Fourier transform were 
obtained by linearly Interpolating the defining points In the magnitude and 
phase of the DFT of the Input and output sequences. The magnitude shown In 
figure 5 Is presented on a decibel scale, normalized with respect to the 
magnitude of the largest component In the Fourier transform. The magnitudes 
of the largest components In the Fourier transforms of the Input and output 
were 2.3xl0~^ and 9.3xl0~7 V- sec, corresponding to frequencies of 0.4 and 
1.4 MHz, respectively. When the magnitude was small, both the magnitude and 
the phase were erratic because of sensitivity to noise contamination (fig. 5). 

Letting the 1024-point DFT of the unit-sample response h(n) of the 
discrete-time system be defined as H(k), we can then use equation (4-4b) In 
table 4 to relate the DFT of the Input X(k) and the DFT of the output Y(k) 
as 

Y(k) = H ( k ) X(k) (3) 

where 0 < k < 1023 for this case. Then H(k) for the system can be evaluated 
by rearranging equation (3) as 

H(k) = Y(k)/X(k) (4) 

The Inverse DFT of H(k) Is the desired unit-sample response h(n). The 
procedure for obtaining h(n) from known x(n) and y(n) Is called 
deconvolution. 
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However, results from attempts to perform the direct division In equation 
(4), using X(k) and Y(k), as shown In figure 5, were unsatisfactory because 
of noise. Thus, the output y ( n ) was filtered digitally before division. 

The digital filter was selected by trial and error. Various digital filters 
were applied, and the resulting H(k) and the unit-sample response h(n) were 
used to evaluate output properties for various Inputs and were compared with 
experimental results. The digital filter selected was a bandpass filter uti- 
lizing time windowing via a 701-point Blackman window (ref. 7). Specifically, 
the digital filter was obtained via the following steps: (1) obtain the unit- 

sample response h(j ( n ) of a desired Ideal bandpass filter of 0.2 to 4.5 MHz, 

(2) shift h(j(n) by 350 points to obtain the unit-sample response h 3 (n), 
and (3) multiply h 3 < n ) by a 701-point Blackman window W](n) to obtain a 
701-point unit-sample response h 4 (n) of the filter. These steps are Illus- 
trated In figure 4 also. The resulting digital filter was a causal, linear- 
phase bandpass filter. Step 1 Identified the desired bandpass filter. Because 
an Ideal filter Is noncausal, as Indicated by Its h d (n), which Is nonzero 
for n less than zero, step 2 shifted h d (n) by delaying It by 350 points. 

A causal filter Is especially Important for real-time signal processing. 

Because the desired digital filter was a finite Impulse response (FIR) 
filter, which by definition should have a unit-sample response of finite dura- 
tion (ref. 7), h 3 (n) was truncated and only numbers between n of zero and 
701, Inclusively, were retained. To minimize the effects of truncation, 
h 3 (n) was multiplied by a 701-point Blackman window w-j(n) In step 3 to 
form h 4 (n). The length of 701 points minimized the effects of truncation 
and yet allowed a 1024-point DFT evaluation. The length of the sequence 
h 4 (n) of the filter affected the DFT evaluation because filtering In the 
frequency domain (l.e., multiplying the responses of the filter and the signal 
In the frequency domain) corresponded to convolution In the time domain. 

Because only the first 80 numbers In the sequence y(n) were nonzero, proper 
evaluation of the linear convolution of the sequence y(n) with the sequence 
h 4 (n), which was nonzero for the first 701 numbers, required at least a 
780-point DFT, according to discussions associated with equation (4-4) In 
table 4. By Increasing the number of points to the nearest Integer power of 2, 
which for this case was 1024, 1024 point sequences and 1024 DFT ' s resulted. 
Thus, sufficient zeroes were appended to any sequence to extend Its length to 
1024 points, unless otherwise Indicated. Also, 1024-point DFT 1 s evaluated via 
a 1024 point FFT were used, unless otherwise Indicated. 

The 1024-point DFT of the unit-sample response of the digital bandpass 
filter designed by time windowing the unit-sample response of an Ideal bandpass 
filter with a Blackman window Is denoted by H 2 (k). Using equation (4 -4b) In 
table 4, the DFT of the output after digital filtering was H 2 (k) Y(k). 

Then equation (4) could be modified as 

H](k) = H 2 (k) Y(k)/X(k) (5) 

where H^(k) was an approximation to H(k), and 0 < k < 1023. 

The unit-sample response h-j ( n ) was obtained from H-|(k) by performing 
an Inverse DFT via an IFFT algorithm. Because h q ( n ) corresponded to h a (t), 
which Is a real function of time, hq ( n ) was a sequence of real numbers and 
corresponded to the real part of the IFFT. Because a shift of 350 points (cor- 
responding to 17.5 ysec) was Introduced In the design of the causal, digital 


320 



bandpass filter, the unit-sample response h-|(n) had to be shifted back by 
350 points to remove the artifact Introduced by the digital filter. This was 
done after the rectangular windowing procedure discussed In the next paragraph. 

The Individual values In the sequence h-j ( n ) for n < 355 and n > 450 
were small, less than l/lOO^h of the maximum value In the sequence. These 
small values were due to noise contamination Introduced during signal recording 
and also by the digital filtering. So, the sequence h-j ( n) was time 
windowed by a rectangular window W 2 (n) defined by 

(1 for 355 < n < 450 

w 2 (n) = < (6) 

(0 otherwise 

Thus, the numbers In h](n) were not affected for 355 < n < 450, but all 
others had been set to zero. The resulting sequence obtained by multiplying 
h-|(n) and W 2 (n) Is denoted as h 2 (n). This time windowing significantly 
smoothed the phase of H-|(k) In the high-frequency region. 

After the rectangular windowing, the sequence h 2 (n) was shifted back 
350 points to remove the shift (l.e., time delay) Introduced by the bandpass 
digital filter discussed earlier. The 1024-point sequence h 2 (n) was shifted 
to obtain h(n) according to circular shifting, which In this case became 

| h 2 (n + 350) for 0 < n < 1023 - 350 

h 2 (n - 1024 + 350) for 1024 - 350 < n < 1023 (7) 

0 otherwise 

The sequence h(n) was the unit-sample response of the test system. The 
numbers In h(n) were zero for n > 100 due to the windowing described In 
equation (6). Circular shifting was used Instead of linear shifting because 
h 2 ( n ) was the result of a few DFT evaluations and, when a DFT evaluation 
was performed, a finite-length sequence was Implied as a period of a periodic 
sequence. Thus, circular shifting procedures corresponding to shifting a 
periodic sequence were required. 

The DFT of the unit-sample response h(n) was H(k). Because h(n) was 
nonzero only for <100, the first 512 points of h(n) were used to compute 
H(k) via a 512-point FFT. Thus, H(k) as obtained was a 512-point sequence; 
H(k) was the frequency response of the discrete-time model of the test system. 

Figure 6(a) shows the Impulse response h a (t) of the test system as 
obtained by linearly Interpolating the defining points In the unit-sample 
response h(n). An Impulse response has units of response per unit of excita- 
tion multiplied by time (ref. 11). Thus, h a (t) In figure 6(a) Is shown with 
units of volts per volt-second. The duration of the unit-sample response was 
approximately 3 ysec. Figure 6(b) shows the magnitude and phase of the fre- 
quency response H a (jft) of the test system as estimated by linearly Inter- 
polating the defining points In the frequency response H(k). A frequency 
response magnitude has units of response per unit of excitation (ref. 11). 

Thus, the magnitude of H a (jS2) has units of volts per volt. The magnitude 
of H a (jQ) Is shown In figure 6(c) on a decibel scale, normalized with 
respect to the magnitude of the largest component In H a (jQ), which Is 
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2.7x1 V/V corresponding to a frequency of 1.7 MHz. The normalized magni- 
tude had a maximum value at 1.7 MHz, and the values of the normalized magnitude 
decreased by 6 dB from the maximum value at 0.6 and 2.3 MHz. The phase of 
H a (j«) was linear from 0.3 to 2.7 MHz. The linear-phase behavior Indi- 
cates that Inputs of frequencies from 0.3 to 2.7 MHz will result In outputs 
delayed In time without distortion. (However, distortion occurred because the 
magnitude of H a (jQ) was not constant from 0.3 to 2.7 MHz.) The Introduced 
time delay corresponding to the linear phase can be computed to be 0.49 ysec. 

As shown by the magnitude plot, the system was Insensitive to frequencies below 
0.6 MHz and above 2.3 MHz; the Information In these frequencies was easily con- 
taminated by noise and also displayed the effects of digital filtering. Thus, 
to summarize, the system behaved as a linear-phase bandpass filter In the 
frequency range 0.6 to 2.3 MHz. 


RESULTS AND DISCUSSION 


The outputs of the test system due to different Inputs can be predicted 
by using the unit-sample response h(n) and frequency response H(k) of the 
discrete-time model of the test system. If the Input Is x(n) and the output 
Is y(n), the Input-output relationship Is given by the convolution In 
equation (1-1) In table 1. For finite-length sequences, equation (1-1) can be 
written In a similar form as equation (2). The convolution Is conveniently 
evaluated via the DFT. If the DFT of the Input Is X(k), the DFT of the output 
Y(k) can be evaluated by using equation (3). So, the output y(n) due to the 
Input x(n) Is the Inverse DFT of Y(k). The frequency response H(k) Is 
shown In figure 6. Figures 7 to 10 show the outputs of the system to various 
Inputs. The time histories of the predicted outputs were obtained by linearly 
Interpolating the defining points In the sequence y(n) as obtained via the 
DFT procedure for 512-point DFT's. 

Figures 7 and 8 show the outputs of the system to tone burst Inputs of 
approximately 6 ysec In duration and 1.5 and 2.5 MHz In center frequency, 
respectively. The agreement between the experimentally measured and the 
predicted outputs was excellent. The arrival times of peaks of Individual 
cycles were predicted to within ±25 nsec In both cases. The measured peak 
amplitudes of all the Individual cycles were correctly predicted to within 
maximum errors of 5 and 9 percent for the 1.5- and 2.5-MHz Inputs, respec- 
tively. Similar agreement was also obtained for center frequencies ranging 
from 0.5 to 3 MHz. 

Figure 9 shows the output of the system to a 7.5-ysec-duratlon rectangular 
pulse Input. The agreement between the experimentally measured and the pre- 
dicted outputs was excellent. The arrival times of peaks of Individual cycles 
were predicted to within ±25 nsec. The measured peak amplitudes of the Indi- 
vidual cycles were correctly predicted to within a maximum error of 5 percent. 
From figure 9, for the approximate time of 5 to 7 ysec, both the measured and 
the predicted outputs were constant In time. The constant values corresponded 
to 0 and 0.046 V for the measured and predicted outputs, respectively. Sim- 
ilar agreement was also obtained for pulse durations ranging down to 0.85 ysec. 
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Figure 10 shows the output of the system to one cycle of 1.75-MHz- 
frequency Input at 2.3 V peak to peak. The arrival times of peaks of Individ- 
ual cycles were predicted to within ±25 nsec. The shapes of the measured and 
predicted outputs agreed well, but the amplitude of the predicted output was 
somewhat larger than that of the measured output. As shown In figure 10, the 
predicted and measured peak-to-peak output amplitudes were 8.7 and 7.2 V, 
respectively. The predicted output was larger than the measured output when 
the Input and output exceeded 1.4 and 5.1 V peak to peak, respectively. This 
nonlinearity Indicates that the linear model used In the prediction becomes 
Inadequate for large-amplitude Inputs and outputs. When the Input and output 
were less than 1.4 and 5.1 V peak to peak, respectively, the measured output 
peak amplitudes were correctly predicted to within a maximum error of 
6 percent. Restrictions on both Input and output amplitudes were Imposed 
because the nonlinearity may be due to the Input or the output or both. 


CONCLUSIONS 

The SWF ultrasonic test system Input-output characteristics were Investi- 
gated by directly coupling the transmitting and receiving transducers face to 
face without a test specimen. Some of the fundamentals of digital signal proc- 
essing were summarized. Input and output signals were digitized by using a 
digital oscilloscope, and the digitized data were processed In a microcomputer 
by using digital signal-processing techniques. 

The continuous-time test system was modeled as a discrete-time, linear, 
shift-invariant system. In estimating the unit-sample response and frequency 
response of the discrete-time system. It was necessary to use digital filtering 
to remove low-amplitude noise, which Interfered with deconvolution calcula- 
tions. A digital bandpass filter constructed with the assistance of a Blackman 
window and a rectangular time window were used. Approximations of the Impulse 
response and the frequency response of the continuous-time test system were 
obtained by linearly Interpolating the defining points of the unit-sample 
response and the frequency response of the discrete-time system. The test 
system behaved as a linear-phase bandpass filter In the frequency range 0.6 to 
2.3 MHz. These frequencies were selected In accordance with the criterion 
that they were 6 dB below the maximum peak of the amplitude of the frequency 
response. 

Furthermore, by using the unit-sample response and the frequency response 
of the discrete-time system and the assumption of a linear, shift-invariant 
system, the output of the system to various Inputs was predicted and the 
results were compared with the corresponding measurements on the system. The 
predicted output was obtained by linearly Interpolating the defining points of 
the discrete-time output sequence. Tone bursts of various center frequencies 
and durations, rectangular pulses of various durations, and angle cycle Inputs 
at 1.75 MHz with various amplitudes were considered. Excellent agreement 
between the predicted and measured results was obtained for tone bursts with 
center frequencies from 0.5 to 3 MHz and rectangular pulses. Excellent agree- 
ment between the predicted and measured results was also obtained for one-cycle 
Inputs of various amplitudes, as long as the Input and output were less than 
5.1 V peak to peak, respectively, at a 40-dB system gain. These results are 
specific for the particular set of components In the test system. 
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The results of this study will be useful In the quantitative analysis of 
the SWF when test specimens are Inserted between the transducers. With known 
test system characteristics, the effects on the output signal due to the test 
specimen having a variety of flaw states can be Isolated. Furthermore, using 
and discussing digital signal-processing methods extensively may facilitate 
developments In CANOE. 
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TABLE 1 Input-Output Properties of Discrete-Time Linear, 

Shift-Invariant and Linear Shift-Invariant Systems 
and Their Corresponding Analogs in Continuous-Time 
Systems . 


Discrete-Time 

System 

System Input-Output 
Properties 

i 

Analog in 
Continuous-Time 
System 

Linear 

If x 1 (n) + y 1 (n) 
and x 2 (n) + y 2 ( n ) > 
then ax^(n) + bx^Cn) 

ay (n) + by 2 (n) 
(Principle of Superposition) 

Linear 

Shift-Invariant 

If x(n) y(n) , 

then x(n - k) -*■ y(n - k) 

Time-Invariant 

Linear 

Shift-Invariant 

If <5 (n) h(n) , 

then x(n) y(n) 

oo 

where y(n) = I x(k)h(n - k) 

k=-°o 

(1-1) 

(Convolution Sum) 

Linear 

Time- Invariant 


* 5(n) is called the 


6<n) 



and h(n) is called 


unit-sample sequence and is given by 
n = 0 


( 1 - 2 ) 

n 4 - 0 

the unit-sample response. 
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TABLE 2 Fourier Representations of Continuous^Time Signals 
and of Sequences. 


Fourier Representation 

Synthesis Equation 

Analysis Equation 

Fourier transform of 
continuous- time 
signal x^(t) 

f 00 

X (t) = X (jn)e jnt dfl 

a ZTT a 

(2-la) 

r 00 *n 

X a (jO) = x a (t)e _jnt dt 

(2-lb) 

Fourier transform of 
sequence x(n) into 
the continuous 
function X(e3 w ) 

x ( n > = W X(e^ t0 )e'* a3n dio 

-it (2-2a) 

X(e jt °) = I x(n)e" jWn 
n=-°° 

(2-2b) 

Discrete Fourier 
series of periodic 
sequence x(n) of 
period N 

N 1 N r 1 . jk(2ir/N)n 

x(n) = tr l X(k)e J 

k=0 

(2-3a) 

sew - T 

n=0 

(2-3b) 

Discrete Fourier 
transform (DFT) 
of finite-length 
sequence x(n) of 
length N into the 
sequence X(k) of 
length N 

, v 1 V Jk(2ir/N)n 

x(n) = m 1 X W e » 

k=0 

for ° <_ n < N-l ; 
x(n) = 0 , otherwise. 

(2-4a) 

„ . N v 1 , . -jk(2n/N)n 

X(k) = l x(n)e 

n=0 , 

for 0 <. k < N-l; 

X(k) = 0 , otherwise. 

(2-Ab) 


TABLE 3 Conditions for the Special Case Where There is 

Exact Correspondence (to within a Scaling Factor) 
of Fourier Transform of a Continuous-Time Signal 
and DFT of a Sequence Derived from Sampling [10]. 


Conditions 

a) 

The continuous-time signal x^(t) must be periodic. 

(2) 

The frequency content of x^(t) must be bandlimited. 

(3) 

The sampling frequency must be greater than twice 
the largest frequency component of x^(t) . 

(4) 

The sampling length must truncate x^(t) with a 


truncation interval exactly equal to one period 
(or integer multiple of one period) of x^(t) . 
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TABLE 4 Convolution Representations of Continuous-Time Signals 
and of Sequences , 


Convolution Representation 

Definition 

Fourier Representation 
of Convolution 

Linear convolution of 
continuous-time signals 
x al (t) and x a2 (t) 

( 00 

x a3 (t) “ x al (T)x a2 (t " T)dT 

(4-la) 

X a3 (jn) “ X al (ja)X a2 ( jn) 
(4-lb) 

Periodic convolution of 
periodic sequences 
xj£n) and x (n) , each 
of period N 

N-l 

x 3 (n) = l x 1 (m)x 2 (n-m) 
m=0 

(4-2a) 

X 3 (k) = X 1 (k)X 2 (k) 

(4-2b) 

Circular convolution of 
finite-length sequences 
x^(n) and x 2 (n) , each of 
length N. The sequences 
are represented as 
periods of periodic 
sequences x^(n) and 
X 2 (n),each of period N. 

N-i ^ 

x 3 (n) = l x 1 (m)x 2 (n-m) 
m=0 , 

for 0 £ n £ N-l; 
x^(n) = 0 , otherwise. 

(4-3a) 

X 3 (k) = ^(k) X 2 (k) 

for 0 £ k £ N-l; 

X 3 (k) = 0 , otherwise. 
This is same as : 

x 3 (k) = x 1 (k)x 2 (k) 

(4-3b) 

Linear convolution via 
circular convolution of 
finite-length sequences 
Xji (n) with length 
and X 2 (n) with length N 2 . 
The length of each 
sequence is increased to 
n 1 +n 2“^ ky appending 
sufficient number of 
zeroes to each sequence. 
The resulting sequences 
are represented as 
periods of periodic 
sequences xi(n) and 
x 2 (n), each of period 
N 1 +N 2 -l 

N x +N 2 -2 

x^(n) = l x 1 (m)x 2 (n-m) 

m=0 

for 0 £ n £ *^+^-2; 
x^Cn) - 0 , otherwise. 

(4-4a) 

x 3 (k) = x x (k)x 2 (k) 

for 0 <_k<_ N 3 +N 2 -2 ; 

X 3 (k) = 0 , otherwise. 

(4-4b) 
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TABLE 5 Significance of Sampling Parameters for Deriving 

a Sequence from Sampling a Continuous-Time Signal. 


Sampling Parameter 

Significance 

Relationship 

Sampling frequency 1/T, 
where T is sampling period 
in sec and 1/T is frequency 
in Hz . 

To avoid distortion 
due to aliasing, the 
sampling frequency 
must be larger than 
twice the highest 
frequency component 
contained in the con- 
tinuous-time signal. 

— < i (5-1) 

2 T ^ ' 

where Q 0 /2 is the 
maximum radian 
frequency of the 
continuous- time 
signal. 

* 

Sampling length NT, 
where N is length of the 
sequence derived from 
sampling and T is 
sampling period in sec. 

To improve DFT fre- 
quency domain 
resolution, the 
sampling length must 
be increased. 

Frequency 
resolution (Hz) 

= 1/ (NT) 

(5-2) 


To avoid distortion 
due to "leakage" 
[23,25] resulting 

from truncating the 
signal, the sampling 
length must be 
sufficient to totally 
contain the signal?* 

NT 2 Time duration 
of signal 

(5-3) 


To utilize the best 
efficiency of FFT 
computation, N should 
be an integer power 
of 2. 

N = 2 V (5-4) 

where v is an integer. 


* If the sequence is to encounter a linear convolution evaluation using 
circular convolution via a FFT algorithm, appended zeroes may be 
required to extend the length of the sequence. 


** If the signal is periodic, the sampling length should be equal to 
the period or an integer multiple of the period. If truncation 
can not be avoided, the truncation effects should be minimized 
by smoothing the signal by proper time windowing [7-10]. 



TABLE 6 


Program Listing for Microcomputer-Based N-Point FFT and 
IFFT Computations in BASIC when N is an Integer Multiple 
of 2. 


30000 

30010 

30020 

30030 

30040 

30050 

30060 

30070 

30080 

30090 

30100 

30110 

30120 

30130 

30140 

30150 

30160 

30170 

30190 

30200 

30210 

30220 

30230 

30240 

30250 

30260 

30270 

30280 

30290 

30300 

30310 

30320 

30330 

30340 

30330 

30360 

30370 

30380 

30390 

30400 

30410 

30420 

30430 

30440 

30430 

30460 

30470 

30480 


rem Subroutine to calculate Fast Fourier Transform (FFT) 

rem or Inverse Fast Fourier Transform (IFFT) 

rem 

rem Data are in fr(real) and fi ( imaginary ) arrays 

rem Calculation is in place, output replaces input 

rem The following parameters must be defined s 

rem nfft# : The number of points (must be power of 2) 

rem ±tr% : 1 for FFT, -1 for IFFT 

rem 

rem Check if nfft# is integer power of 2 

xnfftanfft* 

for i# = 1 to 20 

xnfft=.5*xnfft 5 if xnfft=1 then 30170 
if xnf f t> 1 then 30160 

print M ERROR in subroutine FFT ( nff t% wrong II! ) " : print 

stop 

next i% 

rem Check itr% 

if i tr#= 1 or itrjg— 1 then 30220 

print " ERROR in subroutine FFT ( itr% wrong III ) " : print 
stop 

rem Calculate FFT 
mr$=0 : nn^=nf ft#-1 

for m#-1 to nn% 
nl$=nf f t$ 

nl^=nl$\2 : if mr$+nl$>nn$ then 30260 

mr%=(mr% mod nl%)+nl% : if mr%<=m% then 30300 

tr=fr(m#+l) : fr (m#+1 )=fr (mr#+1 ) s f r (mr^+1 ) = tr 
ti=fi(m#+l) : f i (m#+1 )=f i ( mr#+1 ) : f i ( mr#+1 )= ti 

next m% 
nl#= 1 

if nl56> = nfft$g then 30440 
istep%=nl%+nl% : el=nl# 
for m%=1 to nl% 

arg=3 • 1 41 593*( itr^*( 1 -m^) ) /el : wr=cos(arg) : wi=sin(arg) 
for x%=m% to nff^$ stq.p istep# 

^ 

tr=wr*f r( j^)-wi*f i ( : ti=wr*f i ( ) +wi*f r ( 3 % ) 
fr( j^)=fr( i%)-tr s fi(j#)=fi(i#)-ti 
fr(i#)=fr(i#)+tr : f i ( i%) = f i ( i % ) + ti 
next i% 
next m% 

nl^=istep^ : goto 30320 

if itr$=1 then return 
for i#=1 to nfft# 

fr ( i#)=fr ( ±%)/nff t% : f i ( i % ) = f i ( i % ) / nf f t% 

next i% 

return 
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Signal or Sequence 


Fourier Representation 



(a) 



- 2 ir - ir 0 ir Zir 
- flpT fl 0 T 
2 2 


(iii ) 




x (n) 


• 

1 

• 

i 

Tf«. 


• 

• 

• 

< 

♦ 

• 

4 

L- 0 Sn-I n 


Interpreted 
as T (sec) 


x(n) 


a. 1 


4 


~ i 

♦ ~ T 

♦ 


TtttT 


X(k) 


I li lt 




\n-I k 


DFT X(k) 

Interpreted 
as -j|zp (rad/sec) 


(b) 


Fig. 1 Fourier representations of (a) continuous-time signal and 
sequence derived from sampling, and of (b) finite-duration 
continuous-time signal, sequence derived from sampling and 
the sequence with implied periodicity. 
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Fig. 2 (a) Linear convolution of continuous-time signals and (b) 

circular convolution of sequences derived from sampling the 
signals . 
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Digital Oscilloscope 



Fig. 3 Schematic of experimental components with specimen showing the "stress 

wave factor” test configuration. (All tests in this study were conducted 
with transducers coupled face to face without a test specimen.) 




Unit -Sample Response 
of Bandpass Filter 


700 

Blackman Window 


Fig. 4 Schematic of digital signal processing procedures applied to experimental 
input and output sequences to evaluate unit-sample response and frequency 
response of experimental system. 
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Fig. 5 Time history, magnitude and phase of (a) input signal and 
(b) output signal. 
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Fig. 6 Experimentally determined (a) impulse response and (b) 
frequency response of test system. 
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Fig. 7 Measured and predicted system responses: (a) 1.5 MHz center 

frequency tone burst input and corresponding (b) experimental 
output and (c) predicted output. 
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Fig. 8 Measured and predicted system responses: (a) 2.5 MHz center 

frequency tone burst input and corresponding (b) experimental 
output and (c) predicted output. 
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Fig. 9 Measured and predicted system responses: (a) 7.65 Msec 

duration rectangular pulse input and corresponding (b) 
experimental output and (c) predicted output. 
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Measured and predicted system responses: (a) One 

cycle of 1.75 MHz frequency input at 2.3 Volts 
peak-to-peak and corresponding (b) experimental output 
and (c) predicted output. 







